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ABSTRACT 

Several rare coincidences of scales in standard particle physics are needed 
to explain (i) why neutrinos have mass, (ii) why the negative pressure of the 
cosmological dark energy (DE) coincides with the positive pressure of random 
motion of dark matter (DM) in bright galaxies, (iii) why Dark Matter in galaxies 
seems to have a finite phase-space density, and to follow the Tully-Fisher-Milgrom 
relation of galaxy rotation curves. The old idea of self-interacting DM is given 
a new spin: we propose that the neutrino spin-gravity coupling could lead to a 
cosmic neutrino dark fluid with a an internal energy density varying as function 
of the local acceleration of the neutrino fluid with respect the CMB background. 
We link the Tully-Fisher-Milgrom relation of spiral galaxies (or MOND) with 
the relativistic pressure of the neutrino dark fluid without modifying Einsteinian 
gravity. 

Subject headings: Dark Matter; Cosmology; Gravitation 

It is perfectly conceivable to have a universe with zero cosmological constant and zero 
mass of neutrinos, and in fact it is the most logical construction of the universe to preserve 
the most appealing symmetries of particle physics. Yet symmetry is slightly broken in our 
universe for mysterious reasons: we have a tiny cosmological constant of space-time and a 
small mass of neutrinos. The amplitude of the two are also comparable in terms of the energy 
density, pA ~ (O.OOleV^)"^, and pu = (O.OOOleV^)^ x leV for neutrinos of eV range. Mass- 
varying Neutrinos (MaVaNs) has been proposed, assuming that these two mysteries could 
be related (Fardone et al. 2004, Kaplan et al. 2004), e.g., the mass of the neutrino could 
be due to coupling to a Higgs-like scalar field, whose energy density gives the cosmological 
constant effect. Such a coupling leads to a time- dependent effective mass of the neutrino, 
which is unusual but allowed by physics. 

Given that ordinary neutrinos and sterile neutrinos are hot and cold dark matter can- 
didates respectively, the natural question is to ask if dark matter could also be be unified 
by the same fields. Unification of dark matter and dark energy have been an inspiration for 
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many theories, for example, the Chaphgin gas and the k-essence. Many scalar field theories 
can yield a equation of state parameter w for the dark stress between and -1. The naive 
expectation is that in the absence of interactions, neutrinos have large free-streaming length, 
hence cannot drive structure formation on the small scale. Furthermore the oscillation exper- 
iments favor the neutrino mass around 0.05cV, far from dominating the baryons. However, 
self-interaction in the neutrino sector could modify neutrino's properties. It is known that 
the MaVaNs are unstable to perturbations (Afshordi et al 2006), and in fact behave very 
much like Cold Dark Matter, hence can collapse into neutrino lumps at least on the galaxy 
cluster scale (Mota et al. 2008). 

A harder question is whether the DM effects of the MaVaNs-like sclf-intcracting field 
could also reproduce the TuUy-Fisher relation. Self-interacting dark matter (Spergel & Stein- 
hardt 1999) and its variations have been studied extensively in the past. A fixed cross-section 
per unit mass could help on some scales, but appears to be problematic on other scales. This 
suggests a gradual tapering the cross-section might be helpful. Recently it has been pro- 
posed that the DM and DE effects could be unified by a Dark Fluid described by a single 
vector field (Zhao & Li 2008). The key property of this Dark Fluid is a non-linear pressure 
P depending on the gravitational stress | V$p/(87rG'); this was achieved by taking gradi- 
ent of a time-like unit vector field. Here we propose to link the neutrinos with the Dark 
Fluid of non-linear pressure. The challenge is to give neutrinos the appropriate amount of 
self-interaction on various scales, perhaps making its effective mass increasing on larger and 
larger structure, and so that its property is between that of coUisional fluid and ballistic 
particles free-streaming on geodesies. 

Neutrinos decouple from electrons or matter at Zdec ~ 10^*^ with a thermal momentum 
distribution function 

f{p) = {2nh)-' [eMv) + 1]"' , y = P(^) 0.00018(1 + ^)eV/c (1) 

p{z) 

After decoupling, the expansion of the universe makes this distribution dilate with a charac- 
teristic momentum p{z) or "temperature" ~ 1.95(1 + 2;)Kelvin, or deBroglie wavelength of 
~ 0.1(1 + ^)-^cm. 

It is commonly accepted that observed three flavors of neutrinos and their anti-flavors 
arc the same Majorana particle with mixed states of flavors, helicities and masses. This 
means the particle can start its life from the decoupling of Big Bang as a burst of, e.g., 
ultra-relativisitic neutrinos of pure electron flavor of momentum p{Z(iec) ~ lO^eV, it then is 
decomposed into three mass states of mass < mi < m2 < rris, and changes flavors back and 
forth among electron, muon or tau types. The expansion of the universe cools the neutrinos, 
such that at present day it can have speeds p(0)/mi ~ p(0)/m2 ~ p(0)/m3 ~ 30|^km/s, 
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where p(0) ~ 0.00018eV/c, and we assume the three masses are comparable. The slow- 
neutrinos will be accelerated by the growing localized gravitational potential wells, and likely 
become bound. Observers orbiting inside a galaxy cluster or even inside the solar system 
can sometimes overtake the slow cosmic neutrinos, and identify them instantaneously as a 
flux of anti-neutrinos of reversed helicity. The fraction of bound vs unbound neutrinos or 
anti-neutrinos is a measurement of the potential well x), hence the gradient of the bound 
neutrino number density dxn{t,x) is a measurement of gravity d^^it^x). 

It is interesting if mi = 0, or mi <^ 1712 ~ ma, then a cosmic electron neutrino's wave 
function would be split into two parts, an unbound ultra-relativistic mass eigenstate mi ~ 
plus a bound mixture of neutrinos and anti-neutrino mass eigenstates m2 and m^. So electron 
neutrinos are no longer localized particles, but are a probability distribution over the scale of 
the neutrino horizon, which is only slightly smaller than the light horizon ~ cH^^ by a factor 
~ ^ ~ . J'ff-' While the detailed quantum treatment of the problem is beyond the 

scope of this paper, pressure-like long-range interaction has also been discussed by Pfenniger 
et al. (2007). The bound neutrinos in galaxies will experience a pressure-like interaction, 
mediated by a background of ultra-relativistic neutrino mass eigenstate mi. This might 
create the impression of a modified gravity or modified inertia. 

At a classical level, neutrinos might be treated as a fluid with certain pressure, hence 
the flow of neutrinos is neither on null geodesic as photons nor massive geodesies -^Ua — 
u^dpUa — T'g^^u^u'^ = like the CDM particles. Rather the flow is subject to a pressure 
gradient due to self-interaction, so the flow speed's Lagrangian co-variant derivative -^Ua ~ 
— ^^^daU ~ ln[(l + z){l + 2^{t,x))], where s{n) is the fluid sound speed, which 

is likely a function of neutrino energy density n, whose inhomogeneity tracks the uneven 
gravitational potential ^(t,x). 

What gives mass to active neutrinos is an open question, and it is not even clear whether 
the mass is constant or not. Interaction with Higgs-like field in other sectors or perhaps 
self-interactions in the neutrino sector (perhaps with sterile neutrinos) give an effective mass 
m{z, x) to the chirially left-handed observed neutrinos. This argues that the active neutrinos 
form self-coUisional fluid with less free-streaming. Since the exact mechanism of interactions 
are unknown, one could assume most generally that the effective mass could depend on 
environment. Here we consider models in which the internal energy of neutrinos depend on 
the number density n and its time evolution and gradient Vn. 

Unhke MaVaNs, here we assume that neutrinos have zero mass in the absence of per- 
turbations, and the effects of mass appear only in collapsed systems, like galaxies and the 
solar system. We further assume that the neutrino pressure P is some function of a scalar 
fleld and its temporal and spatial gradient. If the scalar fleld is the temperature T or num- 



-4- 



ber density n of the massless neutrinos, then P = P{T,T ^T,T ^VT). The temperature 
T{x, t) of massless neutrinos varies in essentially the same way (up to a constant factor 0.75) 
as the cosmic photon background, namely, T/Tq ~ {n/nQY^^ ~ ^1 + ^^^f^^ / where 
presently T ~ To = 2 Kelvin on average, and drops with the cosmic time as 1/R{t) and is 
slightly inhomogeneous due to gravitational potential here R(t) is the scale factor 

of the universe, and R = 1 presently. We can replace n/(3n) ~ T/T — ^ R/ R + ^/c^ and 
n~^Vn/3 = T~^VT V$/c^, and assume no preferred spatial directions. The equation of 
state becomes P = P{n,h/n, |V$|), which depends on the density and its rate of change 
and spatial gradient (the local gravity V$|), e.g., P = n'''/(|V$|). The original MaVaNs 
corresponds to a special case / = est, and 7 < 0, so that the neutrinos can be modeled as a 
nearly polytropic fluid with an imaginary sound speed P oc ra''' in an expanding universe. For 
our purpose on galaxies, we neglect cosmic expansion for simplicity, and consider galaxies in 
a static universe with a Minkowski background R{t) = 1, and n ~ hqR'^ = uq 200cm~^, 
which can be treated as a constant. 

To be specific we define {3 to be the sound speed the relativistic neutrino fluid in units 
of the speed of light c, and we adopt a variable sound speed such that 

where we set i? = so that the dimensionless sound speed /? = for the uniform univers^, 
and £"0 is certain fixed low-energy threshold in the neutrino sector respectively, whose physics 
is not specified here. The quantity 1/ [3 carries the meaning of the ratio of gravitational wave 
speed c vs the sound speed c/3 of the neutrino fluid, or the "refraction index" of the neutrino 
fluid, which is assumed to be function of the ratio of the gravitational stress 

vs a characteristic stress hEq = of the neutrino fluid. One can motivate such a correction 
by saying quantum effects must be considered to model interactions of the metric with a 
spin 1/2 neutrino fluid in a box of deBroglie wavelength 1/n ~ I/uq ~ 5mm^; when the 
gravitational energy inside this box is just below or above certain threshold Eq, one might 
expect that anti-particles can disappear or appear (reminiscent of the famous Klein paradox), 
hence can change the energy density of neutrinos and anti-neutrinos, hence the sound speed 
of the neutrino fluid. Interpreting in the picture of the the dielectric analogy of MOND 
(Blanchet & Le Tiec 2008, Famaey, Gentile, Bruneton, Zhao 2007), one would say that the 
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neutrino medium is "polarized" by a variable amount due to a varying gravitational stress. 
The adopted refraction index is such that the perturbations in the stress of the neutrino 
fluid propagate as the speed of light (as in vacuum) in strong gravity, but much slower in 
weak gravity; no propagation or pressure in Minkowski space, where V$ = 0. In the solar 
system, where V$ is much bigger than A, cosmic neutrinos behave relativistically 

We consider a classical Newtonian theory where one minimizes the action S given by 

S = j dtd{xRf + Pg-P), P= (r^Eo)/3^ (4) 

where the relativistic pressure of neutrino fluid is included by a term uEqP'^, and Pg + 
is the Lagrangian of gravitational field and the baryons of density field pb, which couples 
to gravitational potential Here we use the non-covariant formulation of gravity as an 
external field, which is the weak perturbation limit of General Relativity. 

The total Lagrangian now resembles the classical MOND of Bekenstein-Milgrom (1984), 

^ - F(y)) + pb^, (5) 



8ttG 
where 
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y v$ 

Apply the Lagrangian equation by varying the total Lagrangian with respect to $, we 
obtain the MOND-like Poisson equation 

2V [/xV$] = SnGpb, (7) 

where 

p = l-dF/dy=l-{l + ^)-'. (8) 

It is interesting that a relativistic neutrino fluid with a non-trivial pressure can give the 
physics of MOND. To match with MOND acceleration, we can set 

A = 3ao = 3.6 x IQ-^V/sl (9) 

The resulting p function here matches that of Zhao (2007), which are shown to be compatible 
with solar system data and spiral galaxy rotation curves, especially the TuUy-Fisher-Milgrom 
relations. The MOND acceleration scale would corresponds to a mass or energy scale in the 
neutrino sector 

Eo = {3aof/{87rGno) ~ 0.1—^ ~ 2eV. (10) 
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Note that any eV range neutrinos with a normal non-relativistic pressure is insignificant to 
galaxy potential. The non-linear coupling of the neutrino pressure to the metric through 
some unspecified quantum effects is the key. In our picture neutrinos are not localized, this 
can create the effects of a ubiquitous dark matter fluid or modified gravity. 

Let's compare with previous models. The neutrino model here is different from both 
traditional hot dark matter model and from MaVaNs, since we do not require a rest mass for 
neutrinos in the limit of a uniform universe. We could relax the assumption that neutrinos are 
massless in a homogeneous isotropic universe, hence create a model somewhere in between a 
hot dark matter model and pure MOND (Angus et al. 2007). A non-zero mass of neutrinos 
in the absence of perturbations would contribute to the Dark Matter budget in the Hubble 
expansion. Letting the mass change with the neutrino number density would create a dark 
pressure as in MaVaNs. 

There are also important differences of our neutrino model and the classical MOND 
Modified Gravity formulation of a similar effect. Unlike in MOND where is fixed at about 
~ Hq/6, our acceleration scale A in principle varies with redshift, as A oc y/n oc (1 + z)^^"^. 
The fact that the neutrino sector has a mass scale as required by the neutrino oscillation 
experiments, could be the physics of the MOND coincidence puzzle, i.e, why "happens" 
to be comparable to present Hq, i.e., any mass scale in the neutrino sector will set a common 
scale on both dark energy and (galactic) dark matter, so A will track the Hubble expansion. 
A non-zero mass would allow neutrinos to condense into galaxy clusters, the neutrino number 
density n{x,t), hence A{x,t) would be spatially dependent. The MOND effects would be 
made bigger in galaxy clusters, plus the fact that neutrinos add to the total mass of the 
system. This could be a welcome direction to mend the problems of MOND in galaxy 
clusters. This idea is also different from the idea of adding plain massive neutrinos in 
galaxy clusters (Sanders 2003, Angus et al. 2007). However, a rigorous co- variant model of 
the evolution of the neutrino number density involves solving the Boltzmann equation for 
perturbations in a cosmological background, which deserves further study. 

One advantage of our neutrino model here is that it recovers the Tully-Fisher relation, 
hence it is immune from the lasting problems of CDM, i.e., over-prediction of dark matter 
density in galaxy centers and over-prediction of dark matter in dwarf galaxies (Gnedin & 
Zhao 2002 and references therein). While centres of CDM halos can have arbitrary density 
and pressure, the neutrino fluid (cf. eq. [2]) has a limited range of its pressure P = nEoP"^. 
There is a maximum pressure difference of 



iPi - P2I = riEo 
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<nEo (11) 



between the center of a dense galaxy where |V$| ^ A = ^/WnGnE^ and the cosmic average 
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where |V$| ~ 0. This pressure of neutrino appears a key ingredient missing in standard 
Dark Matter models. The present amphtude of this pressure uqEq ~ 0.1 times the critical 
density, hence is comparable to the Dark Energy observed. 

This pressure can be understood as collisional cross-section of self-interacting dark mat- 
ter. Incidentally coUisonal dark matter is also related to the tight correlation of black hole 
mass and velocity dispersion relation of galaxies, which can be understood by a balance 
of hydrostatic pressure (Xu et al. 2007) or the constant refilling of the loss cone by DM 
collisions (Zhao, Haenalt, Rees 2002), or the Bondi accretion of fluid dark matter up the 
mean free path (Ostriker 2000). In short a non-conventional neutrino could relate several 
important puzzles ranging from neutrino oscillation in the solar system to galactic astronomy 
and cosmology. 
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